Mon. Not. R. Astron. Soc. 000, [JHl] (2011) Printed 9 August 2011 



(MN I^TfejX style file v2.2) 



Dense Molecular Gas in Nearby Gas-rich Active Galaxies 



o 

(N 
< 

oo 



o 

6 



Xuejian Jiang*, Junzhi Wang, Qiusheng Gu 

^ Department of Astronomy, Nanjing University, Nanjing 210093, P. R. China 

^ Key Laboratory of Modern Astronomy and Astrophysics, Nanjing University, Ministry of Education, Nanjing 210093, P. R. China 



Accepted 2011 August 7. Received 2011 August 4; in original form 2011 March 24 



ABSTRACT 

We present 3 mm observations obtained with the IRAM 30-m telescope for ten nearby 
gas-rieh active galaxies spanning three orders of magnitude in infrared luminosity. 
Emission lines of HCO+(1-0), HCN(l-O), and C2H(l-0) were simultaneously detected 
in most galaxies of our sample. We also tentatively detected the rare isotopic species, 
Hi^CN, HC^^N and Y{^^CO+ in up to four galaxies (IC694, NGC3690, NGC4258 
and NGC6240). Our estimation shows that the optical depth of HCN is low to mod- 
erate in these galaxies (r ~ 1 — 5). After comparing the intensity ratios of differ- 
ent molecular emission lines (HCN/HCO+, C2H/HCN, and C2H/HCO+) with the 
infrared luminosity (iipj ), we find that in the infrared luminous sample of five com- 
posite (AGN-|-Starburst) galaxies, these ratios vary with Lir: the HCN/HCO'*' ratio 
increases with Lir, which is consistent with previous studies, while there is a slight 
trend in the composite galaxies that both C2H/HCN and C2H/HCO+ ratios decrease 
with increasing Lir. Although our sample is limited, this trend may possibly imply a 
relation between the overall gas density and the infrared luminosity of galaxies. 
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1 INTRODUCTION 

Dense molecular gas plays an essential role in the forma- 
tion and evolution of galaxies. It could supply the mate- 
rial for starbursts (hereafter SB) and for accretion in super 
massive black holes (active galactic nuclei, AGN) in galax- 
ies. It is also critical in understanding galaxy evolutionary 
processes such as galaxy interaction, AGN feedback, etc. 
In return, its properties are strongly affected by these phe- 
nomena and their accompanying radiation fields, such as 
X-ray or UV radiation. Recently, both observational and 
theoretical studies have discussed the relation between the 
environment or activity of host galaxie s and the physical 
or chemical properties of molecular gas ([H cnkcl ct al.l 199ll: 
Krumholz fc Thompsonll2007l: iKrips et ani2008i : iBaan et all 
20081 ). 

Recent studies have developed some molecular gas trac- 
ers, based on millimeter observations, to investigate the in- 
trinsic properties of active galaxies. One advantage of these 
molecular tracers is that their emission lines are less affected 
by dust extinction than optical lines. This is especially im- 
portant for the study of (ultra) luminous infrared galax- 
ies ((U)LIRGs), where substantial amount of dense gas and 
dust are affected by ongoing violent star formation and/or 
a powerful AGN. 

While CO(J=1-0) can trace the general molecular hy- 
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drogen gas content of low to medium density, molecules with 
high dipole moment only trace much denser molecular gas. 
Nowadays, dense gas tracers, such as HCN, HCO^, and CS, 
have been observed and analysed extensively for the study 
of the densest molecular gas content (n(H2) > 10^ cm~^) in 
galax i es (e.g., Geio fc Solomon 2004a| : Gracia-Carpio et al.l 
I2OO8I : iBaan et al.l l2008l : iBavet et al.( |2009| ). A tight Un- 
ear correlation between the infrared luminosity (Lir) and 
HCN(J=l-0) line luminosity (Lp 



HCN{l-0), 



has been found, 



indicating an intimate relatio n between star formation rate 
(SFR) and the dense gas mass (|Gao fc Solomo n 2004b) . This 
correlation has been verified over eight orders of magnitude 
of infrared lumin osity, from the molecular clumps in the 
Galaxy to LIRGs (|Wu et al.ll2005l ). and fr om nearby galax - 
ies to high-z galaxies and QSOs as well (|Gao et al.ll2007l ). 
On the other hand, it has been proposed that the chemi- 
cal and/or physical properties of these molecules vary with 
the types of galactic activity. For instance, HCN/CO(J=l- 
0) and HCN/HCO+(J=1-0) intensity ratios were both found 
to be higher in AGN- dominated systems compared to SB- 
dominated ones (e.g.. lKohnoll2005l : lKrrps et al.ll2008i ). This 
may be due to X-ray dominated regions (XDR) produced 
by AG N which affects the abun dance of HCN and/or 
HCO^ (|Gracia-Carpio et al.l |2006| ). However, it has also 
been demonstrated that there is a wide range of excita- 
tion among different galaxies, implying that only adopting 
HCN(J=l-0) or HCO+(J=1-0) as dense gas tracers could 
be insufficient. Thus there is an urgent need for data of ad- 
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ditional molecular species and transitions, to constrain the 
relati on between dense g as, star formation and AGN activ- 
ities (jGreve et al.ll2009h . Most recently, extensive surveys 
have been focusing on other molecules, such as HC3N, pro- 
viding additional diagnostic tools for the study of abundance 
and excitation prope rties of galaxies (jCostagliola et aD201ll : 
iLindberg et"alj|201ll ). 

In this work we present IRAM 30m single-dish obser- 
vation towards ten nearby active galaxies. HCN(J=l-0), 
HCO'^(J=1-0), and C2H(N=l-0) emission lines were simul- 
taneously detected in most sources. With the simultaneous 
detections and measurements of these molecular lines, we 
compare the intensities of HCN, HCO+, and C2H for vari- 
ous galaxies and discuss the relation between C2H line lu- 
minosity and infrared luminosity of galaxies. 

Throughout this paper, we adopt a flat cosmological 
model with Q.m = 0.27 and Q. a = 0.73, and a Hub ble con- 
stant Hq = 73 km s"^ Mpc"^ (|Bennett et aLlbOllI ). 



2 OBSERVATIONS AND DATA REDUCTION 

We selected ten gas-rich active galaxies (see Table [T]) to ob- 
serve molecular lines at a rest frequency around 87 GHz in 
July 2009 with the IRAM 30m millimeter telescope at Pico 
Veleta, Spain. The EMIR receiver with dual-polarization, 
the WILMA back-end with 2 MHz channel spacing ( this 
corresponds to ~ 6.9 km/s for lines near 87 GHz) and 3.7 
GHz frequency coverage were used for this observation. The 
measurements were carried out using the standard wobbler 
switching mode with beam throws of ±120" and a switch- 
ing frequency of 0.5 Hz. Pointing and focus were checked 
about every 2 hours by measuring nearby QSOs with strong 
millimeter continuum emission. The typical system temper- 
ature at 87 GHz was less than lOOK, including ~ 40K from 
the receiver. We read out one spectrum every 12 minuets, 
which gave an effective on-source integration time of about 
five minutes. 

The data were reduced with the CLASS program of the 
GILDA^ package. A first-order fitting was used to remove 
baselines from all combined spectra. The RMS of each spec- 
trum was also obtained when fitting the baselines. The basic 
physical parameters of our sample galaxies, including origi- 
nal RMS and on-source time, are listed in Table [T] 

For the spectra, we did a preliminary identification of 
molecules refering to frequencies from the NIST database 
Recommended Rest Frequencies for Observed Interstellar 
Molecular Microwave Transitionxs 

The velocity-integrated intensities of molecular lines are 
calculated using 

/ = / Tn^bd?;, (1) 

where Tmb is the main beam brightness temperature, 
and ISy is the velocity range used to integrate the 
intensity. Molecular line intensity in antenna tempera- 
ture is converted to main beam temperature Tmb using 
Tmb = T*Feff/Beff with forward efficiency Feff— 95% 

1 http://iram.fr/IRAMFR/GILDAS/ 

^ http://www.nist.gov/pml/data/micro/index.cfm 



and beam efSciency _Be//=81% at 86 GHz. The Half-Power 
Beam Width (HPBW) in this band was 29". 



3 RESULTS 

Figure [T] shows the entire spectrum of each galaxy with the 
position of certain lines being marked. A more detailed ex- 
amination of the spectrum of C2H(A'' =1 — 0) lines of each 
galaxy can be found in Figure [21 Figure O and 3] show the 
spectra of HCO+( J = 1-0) and HCN(J = 1-0) lines along 
with tentative detections of their isotopes. Table [2] presents 
our measurements of the intensities of molecular lines. 

The results of the prominent molecular emission 
hues HGO+(1-0) and HCN(l-O) for the five (U)LIGRs 
of our sample are consis tent with those observed in 
iGracia-Carpio et all (120081 ). and resuh s of NGC6240 
NGC 4388 and Mrk 231 are consistent with lCostagliola et al.l 
(|201ll ). Meanwhile, we find that some rarely observed species 
toward extragalactic objects are also detected. At > 2a 
level, we tentatively detected: H1^CO+(1-0) in NGC 3690 
and NGC 4258, H"CN(l-0) in IC694 and NGC 6240, and 
HC^^N(l-O) in IC 694. We present the measurements of these 
isotopes in Table U Weak molecular lines are also overplot- 
ted upon the strongest line, with a identical velocity range, 
to be checked if they have a similar velocity distribution. 
Other molecular lines such as HNCO(5-4) and HCO(l-O) 
might be detected in some galaxies. 

There are two groups of hyperfine features of C2H(1- 
0), spaced by ~ 300 km/s and having a line intensity ra- 
tio of 2:1 in the optical thin case under LTE condition 
(jTucker . Kutner fc Thaddeuj|l974l ). Although it is impos- 
sible to distinguish every single component of C2H(l-0) due 
to the line broadening of galactic kinematics, the split hyper- 
fine line features can be readily seen in the spectra of most 
galaxies with prominent C2H detection, and their relative 
intensity ratio also approximates 2 (See Figure [21 where the 
six hyperfine transition frequencies are indicated). However, 
C2H(l-0) of NGC 4258 is an exception, the peak intensities 
of whose split structures are almost equal. 

Generally, C2H(l-0) emission line is the third strongest 
line after HCO+(1-0) and HCN (1-0) within our observation 
window. Determined from the high S/N spectra, its peak 
brightness temperature is typically 1/4 to 2/3 that of HCN 
for our sample galaxies. 

Here we present detailed information of some galaxies 
respectively. 

3.1 NGC 6240 

NGC 6240 is a local LIRG with extreme st arburst activity 
and t he first Binary Active Galactic Nucleus ([Komossa et al.l 
l2003h . The detection of H"CN line was only marginal 
and its spectrum indicates complicated velocity compo- 
nents, probably a refiection of the dynamical structure of 
NGC 6240 (Figure. [4}. Comparing to other galaxies, the 
HCO+ and HCN lines of NGC 6240 show very sharp single- 
peak profile a nd m uch broader line width (> 1000 km s~^). 
iGreve et al.l (|2009f ) reported an extensive molecular line 
data set of NGC 6240, including muhi-transition CO, HCN, 
HCO+, and CS. Our observation of C2H(l-0) lines of 
NGC 6240 is a complement to their work, because C2H(l-0) 
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Table 1. The Basic Properties of Our Sample Galaxies* 



Source 


RA (J2000) Dec 




On-Source time^ 


RMS'^ 




Tvoe'l 




h m s 


o / // 


(km s-i) 


(min) 


(mK) 


(IO^OLq) 




IC 694 


11 28 34.2 


+58 33 48 


3100 


152.5 


1.1 


321 


AGN+SB 


NGC 3690 


11 28 31.0 


+58 33 48 


3180 


111 


1.4 


25I 


AGN+SB 


NGC 6240 


16 52 58.9 


+02 24 04 


7368 


225.5 


1.0 


611 


AGN+SB 


Mrk 231 


12 56 14.2 


+56 52 25 


12500 


133 


1.3 


293I 


AGN+SB 


Mrk 273 


13 44 42.1 


+55 53 13 


11315 


72 


2.0 


125 


AGN+SB 


NGC 4258 


12 18 57.5 


+47 18 14 


500 


165.5 


1.0 


1.2^ 


AGN 


NGC 4388 


12 25 46.7 


+12 39 44 


2560 


171.5 


1.0 


6.1 


AGN 


NGC 4593 


12 39 39.4 


-05 20 39 


2389 


52.5 


1.8 


2.1 


AGN 


NGC 5506 


14 13 14.8 


-03 12 27 


2000 


62 


1.7 


2.9 


AGN 


NGC 5548 


14 17 59.5 


+25 08 12 


5120 


62 


1.3 


3.9 


AGN 



*Tliis table lists our sample. 
■^On-Source observation time. 

'Baseline RMS (corrected to T^b ) of the raw spectrum. 

'^Infrared luminosity Lir from 1) iGracia-Carpio et al.l l|2006h 2) iRice et al.l l|l988l ) and others calculated from IRAS 
infrared flux. 

"^Classification of the galaxies. 



Table 2. Molecular Lines measurements* 



Source 


5V 


RMS'^ 


rHCO + 
■'1-0 




-"l-O 






yHCO+ 


yHCN 




(km s~i) 


(mK) 


(K km s-i) 


(K km s-l) 


(K km s-l) 




(km 


s"i) 




IC 694 


26.892 


0.54 


3.32+0.14 


1.94+0.14 


1.29+0.15 


600 


700 


3096 


3102 


NGC 3690 


26.892 


0.72 


2.15+0.16 


1.10+0.16 


1.20+0.18 


440 


550 


3155 


3165 


NGC 6240 


40.338 


0.41 


4.97+0.18 


2.96+0.18 


1.22+0.16 


1200 


900 


7373 


7363 


Mrk 231 


26.892 


0.66 


1.28+0.15 


1.81+0.15 


0.44+0.18 


500 


675 


12510 


12470 


Mrk 273 


53.784 


0.70 


1.02+0.42 


1.28+0.42 


< 0.25 


800 


445 


11290 


11360 


NGC 4258 


26.892 


0.49 


2.39+0.12 


1.34+0.12 


0.78+0.14 


550 


700 


450 


452 


NGC 4388 


53.784 


0.36 


0.61+0.15 


0.26+0.15 


0.35+0.20 


380 


700 


2534 


2569 


NGC 4593 


53.784 


0.63 


0.19+0.22 


0.34+0.22 


< 0.23 


275 


400 


2394 


2456 


NGC 5506 


53.784 


0.61 


0.44+0.21 


0.20+0.21 


< 0.22 


280 


390 


1972 


1872 


NGC 5548 


53.784 


0.45 


0.08+0.09 


< 0.04 


< 0.06 


100 


200 


5100 


5100 



*RMS and molecular line integrated intensity is corrected to Tmb using Tmb = '^a ^ef f / B^f f where F^f f =95% 
and i?e//=81% in 86 GHz during our observation. Upper limits arc < 2 cr. 
"Velocity resolution of the smoothed spectrum. 

'Baseline RMS of the smoothed spectrum at the corresponding vclosity resolution. Errors of molecular lines is la 
calculated from tr = RMS • VAi; ■ SV , 

'^Velocity range used to integrate the line intensities of HCN and HCO"*". 
''Velocity range used to integrate the line intensities of C2H. 



has an intermediate critical density lying between HCN and 
CO and has its own distinctive chemical property. 

3.2 IC694 & NGC 3690 (Arp 299) 

Arp299, a well-known merging system, is composed 
of two galaxies: IC694 (Arp299A) and NGC 3690 
(Arp299B+C, 45" to the west of Arp299A). Here 
Arp299C refers to the overlapping interface of the two 
mergi ng galaxies, which is also a starburst HII region 



(e.g., iGeh rz. Sra mek fc WeedmanI 1 19831 : lAalto et al.] 1 19971 : 



ICasoU et al. 1999 ). It is believed that IC 6 94 and NGC 3690 
both host an AGN with strong starburst (|Ballo et al.ll2004l : 
ICarcfa-Marm et al.ll2006l'). although f or IC 694 there is still 
some controversy (jGallais et al.ll2004l ). Previous interferom- 
etry observations have revealed that the CO (1-0), ^"^00(1- 



0), and HCN(l-O) emission p eak mainly at the s e three com- 
pone nts Arp299A, B, and C (|Aalto et al.ll 19971 : ICasoU et al] 
Il999h . 

The spectrum of IC 694 shows plenty of molecular lines. 
Besides the obvious HCO+(1-0), HCN(l-O), and C2H(l-0) 
lines, H"CN(l-0), HC^'5N(1-0), HCO(l-O), and HNCO(4-3) 
lines were also tentatively detected at > 2a level (see Figure 
[T]), but H"CO+(1-0) was not detected. 

Of our observation towards NGC 3690, the 29" beam 
only yielded to us a spectrum of Arp299B and Arp299C 
together. The spectrum shows complicated velocity compo- 
nents. No signal of H^''CN line appears, and the velocity 
components of HC^^N are coincident with those of HCN, 
even if the spectrum was less smoothed(see Figure|3)). How- 
ever, considering the low S/N of these two isotopes and the 
high N/^^N abundance ratio (~ 200 - 400), these features 
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could be artificial. The H^'^CO+ line might be blended with 
the hyperfine lines of HCO, as its frequency lies between the 
four hyperfine structures of HCO. We find that the strongest 
HCO+(1-0), HCN(l-O), and C2H(l-0) line components of 
Arp299B+C aU peak at V - 3180 km s"\ which is corre- 
sponding to the starburst region of Arp 299C. We also find 
that the HCO+ and C2H lines in the spectrum of NGC 3690 
have higher peak temperatures (T|j, ~ 12 mK for HCO"'' 
and ~ 4.0 mK for C2H, respectively) than the correspond- 
ing hues in the spectrum of IC 694 (T|j ~ 9 mK for HCO+ 
and ^ 3.5 mK for C2H), but both lines of IC 694 have higher 
integrated intensities. 

Comparing with previous studies of Arp 299, our mea- 
surement of the flux of HCN line for IC 694 (9.8 Jy km s"^) is 
slightly lower than that from lCasoli et af] l) 19991 ) (13.8 Jy km 
s~^), and our measure ment for NGC 3690 (5.5 Jy km s""*") is 
also lower than that of ISolomon. Downes fc RadfordI (| 19921 ) 
(8 Jy km s~^) (We use a 5.9 Jy/K conversion factor at this 
frequency during our observation). These molecules have 
been detected in Sgr B2 giant molec ular cloud and the star- 
forming KL region of Orion cloud (|Turneij[l989l ). In their 
work the peak temperature ratio (T|j(HCN)/Tjj(H^^CN)) 
is ~ 4 for Sgr B2 and ~ 10 for Orion KL, while our results 
show that this ratio is ~ 4 for IC694. And the HCN/C2II 
ratio is ~ 3.6 in Sgr B2, ^ 6.7 in Orion KL and ~ 1.8 in 
IC 694, respectively. 



3.3 Mrk231 & Mrk 273 

Mrk231 & Mrk 273 are local ULIRGs. It is suggested 
that ~ 1/3 the bolometric luminosity of Mrk 231 is con- 
tributed by AGN, and ^ 2/3 by starburst i n the molec- 
ular ring or disk (|Downes fc Solomd 1 19981 : iFarrah et al.l 
I2OO3I ). With respect to Mrk 273, radio observations have 
revealed a n extreme compact starburst toward its north - 
ern region (jPownes fc Solomc|[l998l : ICariUi fc Ta^l2000t ) . 
The HCN/HCO+ intensity ratios of Mrk 231 and Mrk 273 
are larger than unity, whereas HCN/ HCO^ < 1 in other 
galaxi es. Our results are consistent with lCracia-Carpio et all 
(|2008h . and the intensity rati os of HCN/HCO+ an d 
C2H/HCN are aU consistent with lCostagliola et all (|201ll ). 
At > la level, we might tentatively detect HNCO(4-3) in 
Mrk 231. 



3.4 NGC 4258 

NGC 4258 (M106) is an Sab type LINER/Seyfert galaxy. 
There is a n uclear jet related with the anomalous arm 
in NGC 4258 (|van der Kruit. Oort fc Mathewsonlll972l ). Be- 
sides HCO'^ and HCN, the emission of C2H is also strong. 
Other species such as H"CO+ and HCO, were also tenta- 
tively detected. It is intriguing that the proflle of C2H lines 
of NGC 4258 is different from other galaxies. For NGC 4258, 
the emission of the hyperfine components differs only slightly 
in their peak intensities, whereas for other galaxies the emis- 
sion is evidently dominated by the strongest component 
(3/2-1/2 F=2-l). This may imply that the nuclear region 
of NGC 4258 covered by our beam should be optical thick, 
because a optical-thin case requires a 2:1 intensity ration be- 
tween the two hyperfine components, as mentioned aboved. 
But it should be noted that NGC 4258 is the nearest object 



in our sample and our observation only covered the very 
central region (~ Ifcpc) of this galaxy, while the spectra of 
more distant galaxies such as Mrk 231 and Mrk 273 contain 
the emission from entire galaxy. Such difference should be 
considered when we compare NGC 4258 with those galaxies 
of different distance, as the central part of a galaxy is usu- 
ally resided with more dense gas, more active star formation 
and has higher optical depth. 



4 DISCUSSION 
4.1 C2H in Galaxies 

C2H, a representative hydrocarbon molecule, 
was initially detected in some Galactic sources 
iTucker . Kutner fc Thaddeul 1 19741 '). Its relative abun- 
dance in interstellar medium was found to be the same 
order as that of HCN a nd HCO"''. Its dipole moment (0.8-D, 
I Wilson fc GreenI 1 19771 ) is between such high density gas 
tracers and CO. But C2H(4-3) observations toward a few 
Galactic molecular clouds showed that this molecule arises 
from relatively den se gas area, withn(H2) ~ 10* -lO'^ cm 
(|Watt et al.l 1 19881 '). Observations also showe d that C2H 
is widespread over the inner Galactic plane (|Lisztl [T995l l . 
an d was detected in almost all types of molecular clouds 
(Wootten et all Il980l ). The rotational transitions up to 
iV = 11 - 10 (~ 1 THz) of C2H, h ave been me a ,sured 
JCottlieb. Gottlieb fc ThaddeusI Il983l : ISastrv et all Il98ll : 
iMiiller. Klaus fc WinnewisseJ I2OO0I '). The abundance of 
C2H is high in the beginning of massive star formation, and 
then declines in the central region as the hot core evolves. 
Therefore it was suggested to be a go od tracer of th e initia l 
phase of massive sta r form ation (|Beuther et al.l |2008| ) . 
However, IWalsh et all (|2010l '> argued that in star forming 
region, C2H can also trace gas not directly associated with 
active star formation. 

The formation of C2H is mainly determined by 
the formation and evolution of ionized hydrocarbons, 
C+ ^ CH+ (n=2,3,4,5). CH+ transforms into CH4, 
then C2H^ and C2Ht. C2H molecule is formed from the 



dissociated recombination o f C2H j " and C2Hj" (see e.g. . 
Turner . Herbst fc Terzieval I2OOOI : iMeier fc Turiieil l2005l : 



Beuther et"al]|2008l . for details): 



CaHj + ^ C2H -f H, 
C2HJ + ^ C2H + 2H/H2, 



(2) 
(3) 



The sa me chemical precurso r C2H^ is shared by C2H 
and HC3N (|Wootten et al.lll980t ). Other pathways of pro- 
ducing C2H include the photodissociation of acetylene 
(C2H2) and CH2 reacting with carbon atom: 



C + CH2 ^ C2H-f H 



(4) 



On the other hand, C2H molecules are mainly destroyed by 
reacting with oxygen atoms: 



C2H + 0^C0 + CH 



(5) 



For extragalact ic sources, C2H(l-0) emission was first 
detected in M82 (|Henk el. Schilke fc Ma uersberged Il988l ) 
and then in NGC 4945 (Hcnkcl ct al. 1990*). Untfl the lat- 
est molecular line surveys towards nearby galaxies, C2H was 



Dense Molecular Gas in Nearby Gas-rich Active Galaxies 5 



resolved in a third extragalactic source, IC 342, showing dif- 
ferent morphology from other dense molec ular species, such 
as C^*S(2-1), HC3N(10-9) and N2H+(l-0) (|Meier fc Turneil 
I2OO5I ). C2H(2-1) was reported to be the sec ond strongest 
line a fter CS in 2mm observation of NGC 253 (|Martm et alj 
I2OO6I ) . Its abundance was also found to be o f the same order 
of ma gnitude in NGC 1068 as in NGC 253 (|Nakaiima et all 
An extensive survey of C2H(l-0) in external galax- 
ies is obtained by IRA M-30m broadband receiver EMIR 
ICostagliola et al]|201ll '). Our results provide systematic de- 
tections of C2H line in a sample of active galaxies. Although 
C2H is thought to be directly linked to massive star for- 
mation, comparing to our understanding of its behavior in 
Galactic molecular clouds, the global property of C2H in 
galaxies and how it is affected by galaxy activities such as 
AGN and starburst are still ambiguous. Further observa- 
tions of C2H towards active galaxies is necessary for our 
understanding of its properties, which are associated with 
star formation and/or other processes such as XDR or stel- 
lar feedback, especially in NGC 4258 and NGC 4388 without 
active starburst. 



4.2 Isotopic Lines 

As summarized in Table [H we tentatively detected iso- 
topes of HCN and HCO^ in four galaxies for the first time 
(on > 2cr level). The intensity ratios I(HCN)/I(H"CN) 
and I(HCO+)/I(H"CO+) of NGC 3690, NGC 6240, and 
NGC 4258 are all about 40 - 46. In IC694, H1=*C0+ was 
not detected, and I(HCN) /I(H"CN) and I(HCN) /I(HC^^N) 
are only 16 and 20, respectively. Previous observations 
implied that in nearby galaxies, I(^^CO)/I(^'^CO) ratio 
is ~ 10 and ["C]/f^^C1 > 40 jHenkel fc Mauersbergeij 
1 19931 : iHenkel et ail l201oi : iRiguelme et all I2OI0I )." Refer- 
ring to I(HCN)/I(H"CN) ~ (1- 3) x I(i2C0)/I("C0) 
(|Paglione. Jackson fc Ishizukil Il997l l , such result is consis- 
tent with P^CO]/ f"CO] ratio of ^ 19 found in Arp299 
(|Aalto et al.lll99ll ). However, our isotopic line intensity ra - 
tios are significantly larger than ICostagliola et al] (|201ll ). 
this difference might result from that we did not apply Gaus- 
sian fitting to derive the intensities. 

We roughly assume that the filling factor and 
[i2C]/[i^C] abundance ratio for HCN and HCO+ are equal, 
and we also assume that the filling factor and excitation tem- 
peratures of the isotopes of each molecule are equal. Then 
the the optical depth of each isotope is estimated using: 



Table 5. Derived optical depth* 



R = 



1 - ' 



(6) 



where R is the ratio of measured integrated intensities. 
The ratios I(HCN)/I(Hi3CN), I(HCO+)/I(Hi=*CO+) and 
I(HCN)/I(HC^^N) in four sources are listed in Table H 
f'^C]/[ '^C] ratio was sugges ted to be 40-50 for nearby galax- 
ies by iHenkel et all 1 19931 ). For comparison, here we test 
three [^■^C]/[^''C] abundance ratios of 49, 68, and 89 to derive 
the optical depth of HCN(l-O) and HCO+(1-0). The first 
abundance ratio is from investi gation of nea rby star form- 
ing region (|Wang et al.l2009 ,: ,Chin et al.lll999| ). the second : 



local interstellar medium value ([Milam et al 



third is solar system value (jClavton fc Nittleijl200i l. Various 



2005h and the 





49 


68 


89 


NGC 3690: 








t(HCO+) 


0.42 


1.2 


1.9 




u.uuy 


u.uz 


u.uz 










t(HCO+) 


0.12 


0.8 


1.5 


t(H13cO+) 


0.002 


0.01 


0.02 


NGC 6240: 








r(HCN) 


0.40 


1.2 


1.9 


t(H12CN) 


0.008 


0.02 


0.02 


IC 694: 








t(HCN) 


2.8 


4.1 


5.4 


t(H13cN) 


0.06 


0.06 


0.06 






15n 






100 


270 




IC 694: 








t(HCN) 


5.1 


13.8 




r(HCi5N) 


0.05 


0.05 





•Using various [I2c]/["C] (49, 68 and 89) and ["N]/[15n] 
abundance ratios (100 and 270), we estimated optical depth for 
HCN, H"CN, HCO+, Hi3cO+ in four galaxies. HCN/HC^^n 

optical depth is also derived in IC694 (See text for details). 



iLucas fc Liszt|[l99i : ICasassus et"al]|2005l ). and we adopted 
these representative values for crude estimation. The two 
["N]/[^^N] abun dance ratios of 100 and 270 are adopt ed for 
the same reason (I Chin et al ] |l999l : lLucas fc Liszt|[l99i ) . The 
optical depths of each species derived from solving Equation 
(O are fisted in Table [3 

We find that the optical depths of HCN(l-O) and 
HCO+(1-0) are low to moderate (r < 2) in NGC 3690, 
NGC 4258, and NGC 6240, but high (r ~ 3 - 5) in IC694, 
considering optical-thin case always applies for the isotopic 
emission. Our results are consistent wi th the optical depth s 
found for some other nearby galaxies (|Nguven et al.lll992l ). 
Adopting [^■^C]/[^'^C] = 68, we obtain moderate optical 
depth of HCN or HCO^ (r ~ 1), which seems reasonable for 
these gas-rich galaxies. This may be the case, because when 
we use HCN, H^^CN and HC^^N together to calculate the 



optical depth, adopting [^"^Cl/f^C] 



and [i^N]/rN] 



^CJ/l'^^C] rati os have been propo s ed (mainly from CO ob- 
servations, e.g. IHenkel et al.lll993l : iLanger fc Penziaslll99j : 



100, we yield similar results, r ^ 4.1 and r ~ 5.1, for HCN 
in IC 694. We can speculate that in other galaxies without 
detection of isotopes, the isotopic ratios should be larger, 
thus in these galaxies the optical depth of HCN and/or 
HCO+may be low to moderate. 



4.3 Line ratios — dense gas tracer 

Because our observations only collected the emission from 
the inner region of the galaxies matching the 29" beam size, 
before we try to inspect the relation between the infrared 
luminosity and the emission line ratios, we use images of 
Spitzer MIPS 70pim to determine the infrared luminosity ra- 
tio inside the 29" aperture of each galaxies to match the cor- 
responding region of molecular line emission. Although the 
resolution of MIPS 70/im is lower than that of MIPS 24^m, 
the 24/im emission of a galaxy may be contaminated by hot 
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Table 3. Molecular Line Luminosity and Lines Ratios* 





T 1 

HCO+ 
(lOS K km s-i pc2) 


J 1 

HCN 

(10** K km s-l pc2) 


H 

(10^ K km's-i pc2) 


HCN 


C2H 


C2H 


Sources 


HCO+ 


HCN 


HCO + 


IC 694 


L27±0.05 


0.74±0.05 


0.49±0.06 


0.58±0.05 


0.67±0.09 


0.39±0.05 


NGC 3690 


0.82±0.06 


0.42±0.06 


0.46±0.07 


0.51±0.08 


1.09±0.22 


0.56±0.09 


NGC 6240 


9.95±0.37 


5.91±0.37 


2.44±0.32 


0.59±0.04 


0.41±0.08 


0.25±0.03 


Mrk 231 


7.31±0.88 


10.3±0.88 


2.52±1.02 


1.41±0.21 


0.24±0.10 


0.34±0.15 


Mrk 273 


4.71±1.96 


5.92±1.96 


<1.17 


1.26±0.67 


< 0.20 


< 0.25 


NGC 4258 


0.041±0.002 


0.023±0.002 


0.013±0.002 


0.56±0.06 


0.58±0.11 


0.33±0.06 


NGC 4388 


0.20±0.05 


0.08±0.05 


0.12±0.06 


0.43±0.26 


1.39±1.05 


0.59±0.33 


NGC 4593 


0.07±0.08 


0.12±0.08 


< 0.09 


1.84±2.49 


< 0.68 


< 1.25 


NGC 5506 


0.08±0.04 


0.04±0.04 


< 0.04 


0.47±0.53 


< 1.08 


< 0.50 


NGC 5548 


0.09±0.10 


< 0.05 


< 0.07 


< 0.52 


~1.51 


< 0.78 



* Molecular line luminosity based on L' ■ 



:7r/(41n2)e^b/HCN4(l + ^)" 



l|Gao h Solomonll2004al '). Upper limits are < 2 a. 



Table 4. Molecular isotopes* 



Source 


Hi='CO+ 


HCO+ 


Hi=*CN 


HCN 


HC^^N 


HCN 


H"CO+ 


HlilCN 


HCl^N 


IC 694 
NGC 3690 
NGC 6240 
NGC 4258 


0.053±0.05 
0.052±0.05 


40.2±43 
46.3±48 


0.12 ±0.05 
0.073±0.05 


16.8±9.0 

40.6±32 


0.10±0.05 


20.0±13 



*Molccuar isotopic line emission and line ratios of four galaxies with tentative isotope detections 
2a). 



dust emission heated by obscured AGN, thus it may not be a 
good proxy of SFR in this study. And recent Herschel stud- 
ies suggest that the 70/im luminosity is a go od tracer of the 
global SFR even in hi gh-redshift galaxies (|Maanelli et al.l 
l201ll : lElbaz et al1l201l[ ). For the most distant sources (e.g., 
Mrk 231 and Mrk 273), we skip such procedure since they 
were fully covered by the beam size. 

Ratios of HCN, HCO+, and C2H are listed in Ta- 
ble O These ratios are plotted as a function of beam- 
corrected LiR in Figure [S] Intensity ratios HCN/HCO^ ver- 
sus C2H/HCO+, and HCN/HCO+ versus C2H/HCN are 
also shown in Figure [S] 

Despite our small sample and low S/N of some 
spectra, the HCN/HCO"'' ratios span a wide range in 
LIRGs/ULIRGs, which are also AGN-I-SB composite galax- 
ies. However, in other galaxies with Lir < 10^^ Lq these 
ratios remain n early constant at ~ 0. 5. Such dia gram re- 
sembles that of iGracia-Carpio et al.1 (|2006l . |2008| V There- 
after, HCN has been doubted for its reliability as dense gas 
tracer, as AGN activity may enhance HCN in the XDR en- 
vironment, while keeping HCO^ unaffected, resulting in a 
higher H CN/HCO"*" ratio in galaxi es with prominent AGN 
activity (|Gracia-Carpio et al.ll200^ ). On the other hand, as 
a molecular ion, HCO+ is easily destroyed by recombination 
with electrons, and it is demonstrate d that this may caus e 
HCN/HCO"^ to be higher in LIRGs (|Papadopoulodll2007l ). 
Furthermore, the wide excitation range among galaxies ren- 
ders that simply adopting low J transition of HCN and/or 
HCO^ lines are unreliable for determining the properties of 
galaxies. Moreover, the critical densities (ricrit) of HCN and 
HCO^ (of the same J- level) differ by about one order of 
magnitude (|Evanslll999l '). thus they may trace different re- 



gion in the same galaxy. This should be taken into account 
when one compares these molecules as dense gas mass trac- 
ers and determines the influence of galaxy activity on these 
species. Comprehensive multi-transition multi-species obse r- 
vation is needed for such kind of study (|Greve et al.ll2009h . 

The C2H/HCN and C2H/HCO+ ratios plotted versus 
LiR are also shown in panel b) and c) of Figure [S] Since 
we have obtained the integrated intensities of C2H simul- 
taneously, these three molecules can be used to establish 
reliable different line ratios together to study the gaseous 
properties of the galaxies, avoiding both pointing errors and 
beamsize confusions. Panel b) shows a very intriguing rela- 
tion that C2H/HCN seems to decrease with increasing Lir, 
while in panel c) such trend is looser between C2H/HCO"'" 
and Lir. We note that C2H/HCN and C2H/HCO-*" ratios of 
NGC 4258 resemble those of the most luminous galaxies in 
our sample, although its Lir is the lowest. This may be due 
to the beamsize effect that our observation only obtained 
the emission from the nuclear region of NGC 4258 as it is the 
nearest object, while for other distant galaxies we obtained 
more extended gas emission from both the nuclear region 
and the disk. Hence it may imply that the overall molecular 
gas activity of the IR luminous galaxies is more extreme, 
which is similar to the nuclei of nearby normal galaxies. 

The difference between critical density of these three 
molecules can explain such trend in Figure [5] The critical 
densities {ricrit cc fi^v^) of HCN and HCO^ are much higher 
than that of C2H {^l ~ 2.98D for HCN, fj, ~ 3.92D for HCO"^ 
and fi ~ 0.8D for C2H, respectively), so basically they trace 
regions of different density. First, the relatively high density 
gas component can be traced by HCN and HCO^, and the 
intermediate density component can be traced by C2H, thus 
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the C2H/HCN and C2H/HCO+ ratios will decrease toward 
higher Lm, if galaxies with higher Lm, contain more dense 
gas. Second, considering HCN traces denser gas component 
than HCO^ does, the slope of C2H/HCN versus Lir will 
be steeper than that of C2H/HCO"'" versus Lir, as indicated 
in Figure [S] 

Following the above scenario. Figure [6] shows 
HCN/HCO+ plotted as a function of C2H/HCO+ (panel a) 
and C2H/HCN (panel b). In the previous work, ratios be- 
tween HCN(1-0)/HCO+(1-0) and HCN(1-0)/CO(1-0) were 
used to distinguish AGN- and star burst-dominated galaxies 
l|Kohndl2005l : llmanishi et al.ll2009l l . Here we attempt to re- 
place CO(l-O) with C2H(l-0) in such diagnostic diagram, as 
C2H(l-0) reflects medium dense gas content which may also 
be star-forming gas as well. Due to our small sample size 
as well as the large uncertainties for the pure AGN galax- 
ies, we could not separate the galaxies dominated by AGN 
or starburst clearly based on their location in these plots. 
If this is not just due to the large uncertainties, it may im- 
ply that the molecular properties are not distinct in galaxies 
with different activities, but mainly depend on the overall 
molecular gas density of the galaxy. 

In Figure. [71 we explore the relation between molecular 
line luminosity and star formation rate (SFR) indicated by 
Lir corrected for beam effect. Despite the small size of our 
sample, we find that the data points in the three plots show 
different slopes. The slope of our best linear fit to HCN (or 
HCO^) line luminosity (log L') and log Lir is close to unity 
and consistent to previou s results (|Gao fc Solomonll2004a : 
iGracia-Carpio et al.ll2008l ). But in log LcaH - log Lir plot, 
the slope is superlinear, ~ 1.3. This is consistent with the 
scenario that the correlation between star formation rate 
and gas surface or volume density, depends on the criti- 
cal density of the molecular line chosen to trace t he gas 
IIKrumholz fc Thompsonll2007l : [Bussmann et allbOOSl 'l. Thus 
the different slopes may just reflect the different critical den- 
sities of HCN, HCO+ and C2H. 



5 SUMMARY 

We presented broad-band observations in 3 mm band using 
the IRAM 30m telescope for ten active galaxies spanning 
three orders of magnitude in infrared luminosity. The ob- 
servation had benefit of the improvement of the wide-band 
spectrometer onboard the IRAM 30m telescope. Our main 
results include: 

(i) We obtain systematic C2H(l-0) emission measure- 
ments in addition to HCN(l-O) and HCO+(1-0) with si- 
multaneous observation for a sample of 10 nearby active 
galaxies, half of which also show starburst activity. It is an 
important complementarity of recent multi-species investi- 
gation of galaxies, because the chemical behavior of C2H is 
suggested to be directly related to star forming gas. 

(ii) We obtain tentative detections of H"CN, HC^^N, 
and H1^C0+ in up to four galaxies, NGC3690, NGC4258, 
NGC 6240, and IC 694. We discuss the abundance ratio 
and optical depth of ^^C and ^'^C, and find that the op- 
tical depths of HCN and/or HCO^ are low to moderate 
in NGC 3690, NGC 4258 and NGC 6240 but high in IC694, 
considering all the sources are optically thin for the isotopic 
emission. 



(iii) We find that HCN(1-0)/HCO+(1-0) intensity ratio 
spans a wide range in the ultra-luminous infrared galax- 
ies. On the other hand, C2H(l-0)/HCN(l-0) and C2H(1- 
0)/HCO^(1-0) ratios in (U)LIRGs seem inversely propor- 
tional to Lir. We speculate that these results reflect the 
different critical densities of these three molecules. 

(iv) The luminosity of C2H(l-0) shows a correlation with 
Lir in log-log scale with a slope ~ 1.3, while for HCN(l-O) 
and HCO^(l-O) their luminosities are both nearly propor- 
tional to Lir. 
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Figure 1. Full spectra simultaneously obtained from the IRAIVI 30m telescope. The bandwidth of the spectra obtained is ' 
and the RMS and spectral resolution are listed in Table [T] 
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Figure 2. Spectra of C2H(l-0) lines. Green dot dashed lines mark the 6 hyperfine frequencies of C2H(l-0). The intensity scale is 
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Figure 3. Spectra of HCO^(l-O) (black solid lines) superimposed 
with the scaled-up isotopes spectra, M^^CO^ {red dashed lines) 
within the same velocity range. The intensity scales is in unit 
of K. (See the electronic edition of the Journal for a color version 
of this figure.) 
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Figure 4. Spectra of HCN(l-O) (black solid lines) superimposed 
with the scaled- up isotopes spectra, H^^CN^'red dashed lines) and 
HC^®N('6iiie dotted lines) within the same velocity range. The 
intensity scale is in unit of K. (See the electronic edition of 
the Journal for a color version of this figure.) 
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Figure 5. I(HCN)/I(HCO+), I(C2H)/I(HCN) and 

I(C2H)/I(HCO+) plotted as a funtion of Lir for both AGN+SB 
composite galaxies ('ijeci solid circles) and pure AGN('6iMe solid 
triangle). The ratios are in logarithmic scale. The infrared 
luminosities Ljr are corrected for the beam effect (see text for 
details). 
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Figure 6. left panel: I(HCN)/I(HCO+) vs. I(C2H)/I(HCO+); 
right panel: I(HCN)/I(HCO+ ) versus I(HCN)/I(C2H). 




